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ABSTRACT 


A significant part of radar technology has been concerned 
vith developing techniques vhich permit satisfactory opera- 
tion of radar systems vhen moving targets are in the presence 
of relatively fixed targets (clutter). 

This thesis presents the significant aspects of old, 
improved and modern anticlutter techniques, each of vhich is 
described vith its principal advantages and disadvantages. 

A tabulation is made vhich compares each anticlutter 
technique on the basis of its clutter suppression capability. 

Continuing, an application of digital filters is made by 
designing a fourth-order recursive filter vith Buttervorth 
Characteristics, using a computer program. This filter can 
be used in a range-gated MTI (Moving Target Indication) 
receiver to improve its performance. 

Special consideration was given to the AN/UPS-1 and 
AN/SPS-40B,C radars which are used in the Greek Navy in a 
clutter environment characteristic of the Mediterranean 
Sea. Appropriate modifications to their existing analog 


systems are suggested. 
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I. INTRODUCTION 


The conglomeration of unwanted echoes coming from 
fixed targets is defined as "clutter". These echoes clutter 
the radar display and make the recognition of wanted signals 
difficult. The result degrades the quality of radar 
Tuformation. 

There are several types of clutter such as those caused 
by veather, ground, sea, birds, chaffs, rain, clouds, preci- 
pitation and other meteorological phenomena. In addition, 
second-time-around clutter is often troublesome. 

For these types of clutter, there are a number of reme- 
dies for improving the signal-to-clutter ratio. Consequently 
radar systems are configured using one or more of these 
methods. 

The basic signal parameter which is used for distin- 
guishing fixed and moving targets is the doppler shift in 
frequency of the moving target. This frequency shift is given 


by the relationship 


2V 
Tf = _ (1.1) 


d À 
where ës, is the radial relative velocity component of the 
target and X is the wavelength of the transmitted signal. 
The presence of various frequencies is the basis for 
the use of a filter that passes echoes with a doppler shift 


but rejects clutter which has zero or small doppler shift. 





Over the last tventy years, a concerted effort has been 
made to understand the root causes of clutter, its charac- 
teristics, and to find vays for MTI performance to be 
improved. This effort has resulted in nev signal processing 
techniques, vhich together vith certain antenna arrange- 
ments, minimize clutter problems. 

Analog MTI is one technique vhich has been used in radar 
systems for a number of years by the Greek Navy vith 
limited success. Performance has suffered from the inherent 
difficulty of controlling the delay time and system gains. 
The on-line availability has been generally poor. The radars 
operate well when tuned-up, but experience deteriorating 
performance with time. The effective contribution to the 
overall radar system performance has generally been limited 
due to the inflexibility of the hardware. 

However, many of the disadvantages of the analog MTI 
have been eliminated by digital MTI. This is shown in this 
thesis using computer programs to design an optimum filter. 
No attempt was made to implement the filter with hardware. 
However, the computer results are in a form consistent with 
the format needed for easy digital filter design. 

It is expected that when the entire system is built and 
operated, it will have the inherent flexibility and stability 
required to obtain greatly improved performance against 


clutter. 





Tr OLD TECHNIQUES 


In the design of the original radar systems, the main 
effort vas directed tovard the ability of the radar receiver 
to detect the presence of echoes and extract information from 
them. 

Fundamentally, detection vas limited by the presence of 
noise in the receiver. In addition detection of moving 
targets was inhibited by the presence of clutter frequently 
much stronger than the moving targets themselves. 

Under such circumstances the radar became useless. The 
need to overcome this led to the design of additional analog 
subsystems (techniques), which for that time, offered 
considerable detection improvement. A number of these could 
be mentioned, but those few which are still employed in older 
but existing radar systems will be described. 

Part IL covers the following techniques. 


(1) Frequency Agility 
WEE ee Polcarization 
(iii) Sensitivity Time Control - Logarithmic Receivers 


(iv) Instantaneous Automatic Gain Control 
A brief description of each is given, indicating the 


advantages and disadvantages. 


A. FREQUENCY AGILITY 
Frequency agility is a technique vhich enables one to 


achieve a measure of clutter reduction by changing the 





frequency. Coaxial magnetrons permit a change in the basic 
frequency by mechanical means. Signal enhancement in clutter 
by a frequency change is shown in Figure 1 which is an A-scope 
Presentation of a rooftop radar [Refs. 1,2]. 

Part (a) of Figure 1 [Ref. 1] is with the agility switch 
ON and part (b) OFF. The radar was operated to look at small 
targets over a fairly cluttered ground terrain. In this 
case the bandwidth of agility was about 75 MHz. Clearly a 
distinct improvement is seen and the two hard-targets are 
more easily discerned. 

If there is a doubt about the presence of targets, then 
77:27) 5717777 NO Yİ y switch from"ON to OFF can 
usually help resolve the doubt. 

To obtain frequency agility mechanical adjustments are 
made within the microwave tube. These adjustments may be 
as follows: 

(i) By moving the end wall of the coaxial cavity up 
and down slowly. This movement results in broad- 
band tuning. 

(ii) By varying the ring's diameter of the outer 
periphery of the cavity. This form is independent 
of (i) and causes changes in the tuning sensitivity 
ozukucacavity. 

The second technique gives a tuning speed ranging from 80 


to 200 Hz/second and a tuning bandwith up to 100 MHz. 
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Figure 1. A-scope presentations, (a) with and 
(b) without frequency agility. 
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The advantage of this type of tuning is the improved 
performance as demonstrated by its use in operational systems 
today. 

Hovever, there is a disadvantage, vhen precipitation 
extends over a large volume and completely fills the antenna 
beam. Then, the total energy received from all the preci- 
pitation particles vithin the radar resolution cell can be 
quite large compared vith that received from a conventional 
target and can mask or clutter the radar display. 

In this case the amount of clutter with which the 
target echo must compete may be reduced with narrow beam- 
widths and narrower pulse widths if otherwise acceptable. 

The best method, which can be applied to any radar 
system, of eliminating weather effects, such as heavy 
precipitation, is to operate at lower radar frequencies. 

It may be noted that in weather radar applications, the 
designer has the option of operating at the lower frequen- 


cies which give a better S/C ratio. 


B. CIRCULAR POLARIZATION 

When a radar operates at the higher frequencies and 
the target is of a complex nature like an aircraft, it is 
possible to improve the detection capability in the presence 
65:52:02 770 Dy radiating circularly polarized energy 
(Ref. 3). 

Successful detection is based on the dissimilarity 


between clutter and aircraft echoes. With precipitation 
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and raindrops, the backscattered energy does not depend 

on the polarization of the incident vaves. But in the case 
of the aircraft target the energy is dependent on the 
incident polarization. For example, the radar may use an 
antenna vhich is responsive to rotation of the polarized 
vave in one direction. Since symmetrical reflectors reverse 
the rotation of the vave, that echo eneroy vill not be 
accepted by the receiver. 

Circular polarization is generated by tvo methods as 
outlined in Ref. 3 and as shown in Figure 2(a,b). The first 
part shows a quarter-wave plate which converts linear polari- 
zation to circular by passing the wave through it. The 
second part shows a turnstile junction which is a combination 
of rectangular and circular waveguides properly matched to 
generate circular polarized waves. 

Generally with circular polarization the ratio of target 
to precipitation signal is improved by 8 to 25 dB over that 
obtained with linear polarization. Under normal conditions 
an improvement of 15 to 20 dB is expected. It must be noted 
that in all previous cases the ground reflections are 
considered negligible. 


C. S241:549757771Y TIME CONTROL - LOGARITHMIC 
RECEIVER CHARACTERISTICS 


1. Sensitivfey Time Control (STC) 


One of the earliest and most widespread but less 


sophisticated techniques for attenuating the clutter signal 
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is the gain at close ranges. This is knovn as STC İRefs. 2, 
3, 13]. ihe reduction of clutter from close targets is 
achieved by reducing the likelihood of receiver saturation, 
which leads to a programmable receiver gain as a function of 
time so that the gain will be low for targets at close 

range and high for targets at long range. 

The problem of dynamic range limitation of radar 
receivers was overcome by switching on the STC circuit. 

This simply reduced the gain of the amplifiers, usually at 
the intermediate frequency stage. The gain varies as Bo 
where n can take a value between 2 for weather effects and 
4 for land and sea clutter. 

The resulting improvement appears on the PPI, where 
the STC circuit reduces the excessive brightness at the 
center of the display due to the bunching of the radial 
sweep lines. This anti-clutter circuit is used only where 
needed and should be switched off at other times. Conse- 
quently, the use of this technique depends upon an operator's 
judgment and is more helpful the more experienced the 
Operator. 

In addition to the above, this type of anti-clutter 
technique has the limitations caused by its analog circuitry 
and IF stage processing, which limit the radar's clutter 
reduction and target detection capabilities. 

Similar to this is the fast-time constant (FTC), 


another analog circuit. This is normally a high-pass 
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filter (RC), as outlined in Ref. 3. Its use also depends 
onthe Judgement of the radar operator. It is commonly 
used in conjunction with STC and is effective on all 
targets including clutter at all ranges of the radar. 

7 TOgunithmiıC Receiver Characteristics 

The receiver gain, vhen clutter echoes are present, 
can be reduced to prevent relatively large clutter signals 
from saturating the receiver by using a logarithmic 
amplifier followed by a differentiating circuit or a 
pulse-length discriminator [Refs. 2, 3]. 

This technique is used to reduce the effective 
dynamic range of clutter or jamming signals at the threshold 
detector or display. Since clutter varies in intensity with 
range, the use of a log-receiver reduces this and provides 
a clutter-signal amplitude independent of range at the 
output. 

The transfer characteristic of a log-amplifier has 


the form 


“Aa a intl + Dx) (2.1) 


where a, b are constants of the amplifier, x is the instan- 

taneous amplitude of the input and y is the output [Ref. 2]. 
The purpose of the differentiator circuit, which 

follovs the log-amplifier, is to reduce the mean value of 


the noise return or clutter. This is successful only vhen 
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the clutter has a uniform distribution in range for a period 
greater than that of the FTC circuit. 

At this point, techniques l and 2 are compared. 
Sensitivity Time Control (STC) suppresses clutter echoes, 
which enter via the sidelobes much better. STC automatically 
turns down the gain at close ranges, reducing in that way 
sidelobe signals. The log-amplifier does not suppress side- 
lobe clutter which masks the radar display. Both techniques 
have the disadvantage of requiring perfect adjustment and 
skilled operators. 

There are a few applications where these two tech- 
niques may be employed together. Marine radars use log- 
receivers exclusively because they permit strong and weak 
echoes (returns of sea, buoys and boats) to appear on the 
radar scope. 

The block diagram of Figure 3 corresponds to a log- 
amp followed by a pulse-length discriminator (PLD). This 
form of receiver is an alternative solution having the same 
effect on clutter and noise as the log-FTC circuitry, 
rejecting, as well, forms of land clutter and pulsed 
interference. 

The bandpass filter plays the role of the matched 
filter of the conventional linear receiver, and passes 
pulses within fixed frequency limits which are equal to the 
pulse period of the delay line. 

The disadvantage of this technique is the same as 


the log-FTC, in that loss of detectability results. These 


T; 
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losses reach the level of 4-8 dB for single-hit detection 


and 2-4 dB for u-Pulse postdetection integrations. 


D. INSTANTANEOUS AUTOMATIC GAIN CONTROL (IAGC) 

The IAGC technique is a fast-acting form of AGC to 
prevent the saturation of IF amplifiers. It is used in 
most radar receivers and acts in a similar manner to the 
automatic volume control (AVC) in a radio receiver. 

A block diagram is shown in Figure 4, with greater 
detail given in Ref. 3. The control of the gain of the IF 
amplifier is obtained by using negative feedback. Its 
filter has a narrow bandwidth to minimize the effect of the 
extended clutter echoes. 

There is a limitation on the IAGC loop, which can only 
control one or two stages. For this reason instabilities 
sometimes appear in a system with more stages. 

A second technique similar to IAGC is that known as 
detector balanced bias (DDB), which is shown in the diagram 
of Figure 5. DDB differs from IAGC in three basic respects: 

(i) Does not employ a feedback loop. 

(11) Bias is applied to a diode detector instead 

Of the IF amplifier. 

(amy) Employs a delay line in the absence of the 

feedback loop. 

While some of these techniques used singularly or in 
combination are beneficial in reducing the effects of some 


types of clutter, they are not sufficient in modern 


Eo 
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Figure 4. Instantaneous Automatic Gain Control (IAGC) 
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Figure 5. Detector balanced bias (DDB) 
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sophisticated radar systems. More importantiy they are not 
effective in reducing the types of clutter that occur in 


the Mediterranean environment which is the concern of this 


thesis. 
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11: DBOVED TECHNIQUES 


In Norld Nar II a particular class of radar processors 
made its appearance. This type of processor extracted from 
the echo the doppler shift in frequency caused by a moving 
target. Originally, such radar systems employed CV trans- 
mission. The never systems were able to discriminate 
between the echoes of moving aná fixed targets. 

At that time the problem of clutter suppression was 
approached from many directions; "clutter" meaning the 
return echoes from fixed targets in the vicinity of the 
target of interest. Sea, ground, weather effects, angels, 
birds, rain, chaff, trees and vegetation movement are common 
sources of clutter. 

The analog techniques originally developed usually 
employed delay line cancelers or analog filter banks (matched). 
These techniques are still used in large numbers in contem- 
porary applications. Delay line cancelers are used in the 
UPS-1 and AN/SPS-40B,C radars on ships of the Greek Navy. 

Thus, this thesis continues to compare the performance 
of other techniques, both existing and modern, in an effort 
to find a satisfactory clutter rejection filter that is not 
unduly complex and can be rather easily implemented. 

Part III describes the following representative tech- 
niques with their respective advantages and disadvantages. 


er) Analog MTI and Pulse-Doppler Systens. 
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(a) Dual Beam Antenna. 
7111) Log - CFAR. 


(iv) Video Detector. 


A. ANALOG MIT AND PULSE-DOPPLER SYSTEMS 

When doppler information is extracted from a pulsed 
radar return, it is known as an MTI radar. The MTI stands 
for "moving-target indication". It may also be called 
Pulse-doppler radar. In practice, a distinction is some- 
times made between the MTI radar and the pulse-doppler radar, 
although they are both based on the same principle. 

MTI usually refers to a radar in which the doppler 
frequency measurement is ambiguous but the range measurement 
is unambiguous. Another characteristic feature of the MTI 
radar is the delay-line canceler used to detect the doppler 
IEeguency shift. 

In the pulse-doppler radar the doppler measurement is 
usually unambiguous and the range may or may not be ambiguous. 
Ambiguous range means that multiple-time-around echoes are 
possible, while ambiguous doppler implies that "blind speeds" 
fall within the range of expected target speeds. 

MTI radar has the ability of extracting the moving-target 
echo from the clutter echo even when the clutter is 20 to 
30 dB greater than the (a ss pss echo. Some pulse 
doppler radars can detect moving targets even when the 


57706: ceci is 70 to 90 dB greater than the target echo. 
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The distinction betveen the tvo radars arose histori- 
cally, and is commonly still made. In many instances the 
difference betveen them is only a matter of nomenclature. 

Figure 6 shovs the simplest MTI processor, the single 
delay-line canceler, vhich extracts the doppler information 
by subtracting tvo successive echoes. The delay-line 
canceler plays the role of a filter, eliminating the d-c 
component of fixed targets and permitting the a-c components 
Of moving targets to pass. 

There are several types of MTI processors. A distinc- 
tion betveen them is the type of information processed in 
the returned signals, that is, vhether the phase, the 
amplitude, or both phase and amplitude are processed. These 
types are described vith their block diagrams in Refs. 2 and 
3. It is beyond the stated intention of this thesis to 
describe all types of MTI, but some are listed here for 
reasons of interest. 

(i) Phase-Processing MTI 

(11) Amplitude-Processing MTI 

(iii)  Vector-Processing 

(iv) Clutter-Locking MTI 

(v) Multiple-Canceler Systems 

These processors do not satisfactorily handle the well 
known problems of clutter precipitation. This problem 
becomes more important when the effects of wind shear and 


average clutter velocity are taken into account. For example, 
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the clutter performance of the aforementioned three radars 
is poor under these conditions. This has been verified by 
the experience of the Greek Navy in the Mediterranean 
environment. 

To get a clearer idea of this limitation on MTI per- 
formance, the clutter improvement factor I was plotted vs. 
mange. Figure 7 shows the plot. The equation of I is 


as follows 


| r 

EM ——— — ———”- CRD 
22 777 

8 (cr 5515075”) 


and the Continuous System Modeling Program 1 vas used vith 


some typical values given: 


c variance of wind turbulence, varying from 
.-/ to l.l meter/sec 
À = .30 m, transmitted wavelength 
f, = 600 Hz, PRF 
k = 5.7 meter/ (sec) (km), is the component of wind 
иә иле аси ли the direction of the 
radar beam 
R = range in km 
$ = 1.5°, vertical half-power beamwidth-two-way 


path in radians. 
Figure 7 shows that the attenuation of clutter is quite 


good for ranges up to 18 km where, for o and 


turb Ы 
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Figure 7. Clutter improvement factor vs. Range 





1.1 meter/sec a difference of 1.5 dB is measured. As range 
increases, the clutter attenuation decreases rapidly as 


indicated. 


B. DUAL BEAM ANTENNA 

The MTI visibility is highly degraded at short ranges 
due to nearby graound clutter, sea clutter and angel echoes. 
To eliminate these problems a dual beam antenna system has 
been designed. 

This system is shown in Figure 8. Using the low and 
high beam horn at the appropriate time, as outlined in 
Ref. 13, the clutter is reduced at the receiver input. The 
change in the use of horns occurs by switching. The two 
horn noses have a difference of about 4 degrees in elevation. 
The angle is variable. 

The degree of clutter suppression and target visibility 
can be determined by the beam angle and switching distance. 
The coverage is varied by varying the angle difference in 
the nose directions. Consequently, the switching distance 
1s changed and this determines the degree of clutter 
suppression. 

Figure 9 [Ref. 13] compares the S/C ratio of dual beam 
antennae with single beam antennae. After determining 
target/clutter ratios by a computer simulation program, 
conclusions from these two techniques are obtained. This 
comparison shows that a better improvement factor is 


obtained using a dual beam antenna system. 
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Figure 8. Dual beam antenna system block diagram 
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Figure 9. Comparison in S/C ratios betveen dual beam 
and single beam antenna systems. 
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Types of clutter, vhich are taken into account here 
are ground, sea and angel echoes. In all cases an improve- 
ment of 15 to 20 dB is obtained using dual beam antennas. 
Detection, however, becomes difficult when the range is less 
than 5 nautical miles. 

These advantages of the dual beam antenna are the reasons 
which make it useful in systems such as ARSR (Air Route 


Surveillance Radar). 


C. LOG-CONSTANT FALSE ALARM RATE 

Additionally, a special technique for ARSR systems is 
that of Log-CFAR, which has been widely used to eliminate 
weather clutter. The reason this technique has been 
employed is that neither MTI nor dual beam antenna systems 
have the ability to suppress weather clutter which has a 
radial velocity distribution at high angles of elevation. 

Previously mentioned methods obtain some improvement 
by improving the signal-to-clutter ratio. But in the case 
of clutter, whose signal strength is well above the threshold 
or noise level, the output of the receiver is saturated. 

An additional suppression of clutter is needed and is 
obtained by using the Log-CFAR method. It keeps a constant 
output at the receiver suppressing the clutter to noise 
level. 

In terms of the block diagram this technique is shown 


in Figure 10 and its operation is outlined in Ref. 13. A 
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step by step explanation is given in terms of input-output 
vaveforms according to the diagram of Figure 10. 


The output v of log-amplifier is given by 


w = a ln(lzdbv) (SIR) 


where a,b are constants and v is the input [II-C-2]. The 
subtractor has an output y, which is given by 


y = w-x = aln (4) +” (3.9) 


where x is the average detector value given by 


x = = J w, (RIO) 


and where y is the Euler's constant, o the variance of the 
input signal v and N, the sample number. The output signal 


of the anti-log converter is given by 


z = c exp (d-y) = — exp (Z) (3.11) 


where c,d are constants and a.d = l. 

Sample number, N, plays an important role, since it 
determines the processing loss of the CFAR circuit and 
consequently the clutter suppression. Reference 13 shovs 


that the cutoff frequency of the rejection band is given by 
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JE (3.12) 
where t is the pulse width. 
The Log-CFAR process has two disadvantages which are 
listed below. 
(i) Sample number N must have a value of about 8 
to keep the processing loss between l and 2 dB. 
(ii) Mean value of clutter component must have a 


frequency less than 1/Nt to be eliminated. 


E VIDEO DETECTOR 

This application is specially useful in non-coherent 
airborne radars operating either over land or in the proximity 
of a land mass, where very intense ground clutter patches or 
individual scattering objects exist. 

Such ground clutter may have radar cross-sections on the 
order of one to ten million square meters. Thus, even with 
the best designed antennas on airborne platforms, this 
clutter may be strong enough to be detectable through the 
antenna sidelobes. 

These signals in general will have doppler velocities 
within the acceptance band of the MTI filter or filter bank 
and will result in a high rate of false alarms. The side- 
lobe clutter problem is a very serious one but the use of 
a Square-law detector provides significantly improved clutter 


тедес2топ. 
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For main lobe clutter vhich is assumed to be Gaussian 
in both amplitude and frequency, the MTI improvement factor 
is calculated for different spectral vidths using linear 
and square lav detector characteristics. 

The block diagram of the system is shovn in Figure 11 
vhere a magnetron is used as the transmitter. The echo 
signal is converted to IF and then amplified. Detection is 
realized by the so called Video Detector. After MTI 
filtering, digital signal processing follows before the 
video presentation on a display. 

The signal x(t) has a general spectrum in the IF shown 
in Figure 12. Assuming that the clutter component is much 
larger than the other components, the video spectrum is 
shown in Figure 13. After the video detector this spectrum 
is approximately the auto-convolution of the input spectrum. 
Convolutions not containing C, vhich stands for clutter, are 
neglected from the spectrum. Component C*C is undesirable 
clutter and is suppressed by a digital high-pass filter. 

The system described above, using the square law 
detector provides a wider region of good clutter suppression. 
It is clearly superior in comparison with other improved 
clutter suppression techniques, but the dynamic range of the 
output signal is, of course, larger for the square law 
detector which is a disadvantage in several applications. 

It can be noted that at ranges where there is no clutter 
the target is not detectable after the MTI filter. For 


target detection, the filter must be bypassed at such ranges. 
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Clutter beyond the unambiguous range, after convolution vith 
first time around echoes, results in vhite noise. This is 


an additional system problem, 
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IV. MODERN TECHNIQUES 


This section is prepared to demonstrate several 
"modern techniques" vhich are employed in radar system 
technology. The word "modern" covers the developments of 
the past 5 years and includes future trends in signal wave- 
form design and processing technology. It describes three 
techniques, which are listed below. 

(1) Pulse Doppler Processing 

(11) Moving Target Detector 

(iii) el ication of PET Algorithm. 

Finally, a tabulation is made to compare the amount of 
clutter suppression achieved by each technique in the three 


categories on different types of clutter. 


E SRBSE-DOPPLER PROCESSING 

Here a technique is presented which is especially 
applicable to high-power, pulse-doppler radars, achieving 
good clutter discrimination. The basis of this method is 
an adaptable filtering and integration system which processes 
the signal in a way that the processin loss is minimized. 
Besides, good range resolution and velocity resolution can 
be obtained. The integration is made adapative to the 
changes in the variable PRF by using commutated-capacitor 


filters. 


40 





The clutter filter system has a block diagram shown 
in Figure 14. The signal waveform looks like a Gaussian- 
modulated envelope due to its antenna characteristics. A 
Gaussian filter is placed at the beginning, which broadens 
the duration, 1, Of the pulses and consequently reduces 
their peak amplitude. 

This filter provides a narrow spectrum of the signal 
for adjacent channel rejection during the mixing processes. 
Then a Gaussian low pass filter rejects the fundamental 
frequency fo and all harmonics generated by the previous 
process. After the sample and hold circuit a high pass 
filter is connected. The purpose of this filter is to 
attenuate clutter components to below noise level. 

The two signals are then mixed with fe Since these 
are out of phase by 1/2 radians, they are combined to provide 
27:51:000 “ignal from vhich elutter has been removed. 

Some important advantages of this technique are listed 
below [Ret. 16]. 

(1) Clutter suppression at levels of about 80 dB 

down in amplitude. 

(11) Use of capacitors in association with 
adaptable sample and hold circuits, stationary 
clutter is reduced. 

(111) Rapid changes of PRFs during the time on 
target permits the system to respond without 
spreading the clutter or destroying the 


integrity of the signal information. 
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Figure 14. Clutter filtering system 
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(iv) The lov-pass filter reduces system noise 

to a much lover level. 

The disadvantages of the systemare the range and 
velocity ambiguities which occur when it is required to 
detect very fast targets. Also, the greater complexity 
due to three distinct stages of spectrum analysis. 

This study leads to the conclusion that the introduction 
of clutter filters, which are switched in synchronism with 
the pulse repetition frequency sequence, effectively reduce 
the generation of clutter side-bands by an additional 
40-50 dB over a fixed filter system. It is thus possible 
to design a long range Pulse Doppler Radar System with more 
than 80 dB clutter suppression, if the PRFs can be switched 
at a rapid rate. 

The storage elements are capacitors which are switched 
by a hardware program to integrate or discharge signals. 

By this means real-time target assessment is realized. The 
performance of this radar system is consequently significantly 


improved. 


В. MOVING TARGET DETECTOR 

Recently, the Lincoln Laboratory, associated with 
Massachusetts Institute of Technology, developed a signal 
processor named the Fast Digital Processor which performs a 
Digital Fourier Transform on the received signal and presents 


De Pourier coefficients [Refs. 1, 17, 18]. 
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This processor is very fast, since it is a special- 
purpose computer designed to perform signal processing tasks 
such as FFTs. Its algorithm was tested by embodying it in 
a hardwired version of the Moving Target Detector (MTD). 

The MTD has been designed, built and tested at Lincoln 
Laboratory where it is now undergoing further tests. 

The MTD technique uses the cascade of an MTI and pulse 
doppler processor to achieve clutter rejection. This 
combination is an approximation to an optimum clutter 
rejection filter. 

In terms of the block diagram, the MTD is shown in 
Figure 15, where it is seen that the signal to be processed 
arrives at the IF stage from the output of the IF pre- 
amplifier and is fed through a special wide-dynamic-range 
amplifier to the quadrature video detectors. 

The two quadrature video detector outputs are converted 
to 10-bit digital numbers by the A/D converters. The system 
contains an 8000-word input memory and about 900 integrated 
circuits. A disc memory is used as a fine-grained ground 
clutter map. This memory allows the removal of such 
returns by the use of adaptive thresholds in range, doppler 
and angle. 

It is required that a number of scan returns be observed 
before a target is declared present and track is initiated. 
Thus, track initiation is used to help eliminate false 


clutter returns. Multiple’ PRFs are used to detect weak 
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targets that ordinarily vould be masked by rain clutter 
having the same ambiguous doppler frequency. 

Because ground clutter is so spotty in nature, great 
difficulty is experienced in calculating appropriate 
threshold values for detection. To solve this problem a 
digital ground clutter map is implemented vith one vord 
for each range-azimuth cell. All these are stored on the 
magnetic disc memory. Sometimes magnetic-bubble and EBAM 
(Electron-Beam-Addressed Memory) are used, both of which 
have the advantage of being nonvolatile. Also, CCD 
memories are used, in place of disc systems, which are 
Suitable for storing the clutter map. 

MTD is said to have subweather visibility. In fact 
it can see a small commercial type aircraft whose cross- 
section is many dB below the radar cross-section of the 
weather return. This feature was not previously available. 
Now it is used to provide excellent tracking of aircraft in 
rain. 

A second feature of the filter bank approach is the 
facility with which proper thresholds can be established 
taking into consideration the presence of rain. The 
detection threshold is established by summing the detected 
Output in range cells of interest. 

The MTD has the advantage of being able to detect and 
track zero doppler aircraft targets while ignoring bird 
and insect targets. At zero doppler, aircraft targets are 


broad-side and hence have much larger cross sections than 
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do birds or insects. Therefore, a higher threshold can 
be used in a zero doppler filter to eliminate bird echoes 
Bat Stil] detect aircraft targets. 

Tracking is automatically initiated and the tracks are 
continuous despite ground and weather clutter and birds. 
The radar now can be sited freely without consideration to 
g ound clutter limitations. 

Both MTI and pulse doppler processing have the disadvan- 
tage in regard to vulnerability to jamming. For both 
techniques many more than one pulse detection is required 
to see a signal imbedded in clutter. 

Finally, no tuning or adjustment of any kind is needed 
for this combination. It is a reasonably economical 


solution to the radar automation problem. 


C. FFT APPLICATION 

Within the past several years, digital processing of 
radar signals has become quite popular. Recent advances in 
A/D converters, integrated technology, and the emergence of 
the Fast Fourier Algorithm have made practical the digital 
Srltering Of Signals. 

The application of the FFT algorithm starts from the 
first known technique of Convolution by the Discrete 
Fourier Transform. It is shown in Figure 16, where a 


block diagram explains the process. 
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Figure 16. Convolution by discrete Fourier Transform 
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It is simple in concept and consists of For ring rantinpi 
sequence of samples from the received signals. Then a 
Fourier Transform is performed to yield samples in the 
frequency domain, multiplying these by a stored filter 
function and performing an inverse transform, yields the 
time domain response. 

Thus, any filter function can be implemented by storing 
it in memory. So, knowing the clutter spectrum one can, 
in principle, filter the clutter by implementing the 
appropriate function. 

The FFT algorithm is shown in the pipeline of Figure 17. 
This technique is used to digitally realize a bank of 
contiguous filters and thus produce a spectrum analyzer or 
doppler filter bank. 

Advantages of the system include stability, reproduci- 
bility, and flexibility. The amount of time for this 
processing is a parameter of considerable interest because 
real time solutions are commonly required. To reduce time, 
high-speed A/D converters are being developed and employed 
in modern systems. 

Following the discussion of the several techniques 
described in Chapters II, III and IV, a tabulation in 
matrix form comparing filter characteristics with clutter 
characteristics has been prepared and is shown in Table I. 
The classifications of techniques are qualitatively identi- 
fied as Superior, Intermediate and Inferior, in terms of 


the amount of clutter suppression or improvement factor in 
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dBs and also the range at which each process is reasonably 
sensitive where data information is available. 

The purpose of this comparison is to help in making a 
decision about the system to be selected for a give use. 

It is noted that the decision will be baseá on the expected 
requirements from the system the area where it is to be 
installed, if a fixed or mobile site, and on other condi- 
tions which need consideration. 

Each radar system is generally designed to optimize the 
realization of a specific set of measurement or observa- 
tions. Once the design criterion has been achieved, usually 
it becomes difficult to utilize the radar for other purposes 
without extensive and costly modifications. 

The observations and measurements shown in Table I were 
obtained from a number of independent observers. The data 
was taken chiefly at X-band on appropriately configured 


Navy and air-surveillance radars. 
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V. AN APPLICATION OF A DIGITAL FILTER 


In the radars mentioned earlier, use of analog MTI 
systems give a clutter suppression of 13.8 to 24.6 dB for 
rain clutter according to the data of Table I. This suppres- 
sion is not sufficient and the problem becomes greater vhen 
the radar system is operating under Mediterranean environ- 
mental conditions, where the fall and winter period is full 
of heavy rainfalls. 

The experience of the Greek Navy for certain clutter 
due to sea and birds does not concur with the suppression 
values given in Table I. A level of 18 to 25 dB is achieved 
under perfect radiation conditions. 

Since the radars under consideration are in use in the 
Greek Navy, a requirements exists for the improvement of 
their signal processors in order to get an improved clutter 
suppression. 

This improvement can be achieved by the use of a higher 
than second-order doppler filter according to Ref. 12. 
Mens the motivation for the design of a fourth order 
recursive digital filter. 

For this purpose a Butterworth characteristic band- 
pass filter is to be used with an effective bandwidth of 
Мису Лес Certain restrictions are put on the 
filter design as they are defined in the design part of 


this thesis. 


D) 








The analysis and design are done according to Ref. 10. 
The concept of predestorted or prevarped analog transfer 
function is used. To get a one-to-one mapping in the 
z-domain approximately the same as that of the s-domain, 
the well-known bilinear z-transformation is used. 

To detect moving targets, the system uses the doppler 


shift in frequency, which is given by the relation 


7-5 dı) 


The wavelength, A, of the transmitted pulse might also be 
used aS a parameter. Since the design of the above filter 
refers to specific systems, where A is fixed by design 
parameter, this option is not open. These fixed values of 
A for the two radars are calculated by using the relation 


A = 55: 


— 
fo 


where 

c = 3x10° m/sec, velocity of light 

fo - operating frequency of the system given in Hz 
For ə = 1300 MHz (AN/UPS-1) A = .23 m; 
and Lor Ee = 500 MHz (AN/SPS-40BC) A .6 m. 


“ 20046 sor relative velocity of the target. 
Usual types of moving targets, in the area of interest, 


are given in Table II with their cruise and maximum 
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Type Cruise Maximum 
| Helicopters Generally 80-130 180 
| C-47 150 250 
pops 200 300 
ALBATROS 200 350 
Commercial Air Lines 300-400 500 
mu. 250 600 
A-7H 300 600 
| F-S4F 550 600 
Е-5 400-550 700 
F-102 500 800 
F-104 400-550 1000 
F-1CG 500 1000 
F-4 500 1300 
EZELE IL. NELOCITTES OE SELECTED TARGETS 
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velocities in nautical miles per hour. This information 
ies used in the beginning of the filter design in order to 
evaluate the necessary doppler frequencies assuming 

(1) Radar stationary 

(11) Radar in motion with a velocity varying from 

SR ORSO nautical miles per hour, to correct the 
effective bandwidth of the filter. 

Table III contains a summary of the literature [Ref. 2, 
2550 6) Ref. 3, pp. 146, 551-554, Ref. 19, p. 17-9) for 
the relative velocities and standard deviations, Су, of 
clutter spectrum in meters per second. These kinds of 
clutter are common to the Mediterranean environment. 

The conversion of the standard deviation of the clutter 
power spectrum, er to frequency is obtained by applying 


the relation 


g = (5.2) 





Where A is the wavelength of the transmitted pulse in meters. 
This relation is an alternative to equation (1.1). 

Computer programs 2, 3 and 4 were used to plot the 
amplitude square characteristic versus frequency, uw, of 
(1) analog, (2) prewarped, and (3) digital transfer functions. 
Figures 19, 20 and 21 show the characteristics of the designed 


filter. 
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A. DIGITAL FILTER DESIGN 

The second part of this thesis is the design of a 
digital filter which is going to reject clutter echoes in 
ranges indicated by Table III [Refs. Sor SLT Certain 
numbers of doppler filters combined with range gates 
constitute the video part of an MTI radar which is shown 
HEEN 18 in block diagram configuration. 

A good performance criterion of MTI doppler filters is 
the MTI improvement factor. This factor defines the improve- 
ment of Signal to Clutter ratio which can be achieved by a 
filter. In the table of Ref. 12, the improvement factor I 
15 given as a function of the order of the filter. 

Generally the higher the order of the filter the better the 
improvement factor. However, higher order transfer functions 
imply complexity, low reliability and higher cost. 

The goal of the present doppler filter design is to 
derive a transfer function which will provide high clutter 
suppression and at the same time reduce the complexity 
caused by the higher order filter, to a reasonable level. 
For the above reasons, a fourth order digital filter will 
be designed. The improvement factor will be increased and 
the filter complexity will be kept to a minimum 

An Infinite Impulse Response (IIR) band-pass filter of 
the fourth order will be implemented. As is known, the IIR 
filter does not need as many delays as the Finite Impulse 
Response (FIR) filter. Consequently, the number of weights 


is considerably lower. 
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The transfer function for a fourth-order band-pass 


digital filter is given by the following equation 


-1 =”. -3 - 
4- 
H(z) = A(z) . xa o S EE MER (5.3) 
Biz) -1 -2 -3 -4 : 
1 +b,2 +b,2 * 642 +b,z 





This transfer function is completely defined by the 


= 
values of ag! ај, 22, EL PL bi: Do, 3” 


tvo vays to define the above coefficients. The first is 


b Du: There are 

to go directly to the z-domain. According to the frequency 
response characteristics which will be defined later, one 
has to find the exact position of zeros and poles of H(z). 
Since zeros and poles are calculated the values of a's and 
b's occur after a simple multiplication of the factorized 
polynomials A(z) and B(z). 

The second method, to be used here, is to select the 
desired H(s) transfer function. From this is defined the 
transfer function H(z). It is more difficult than the first 
method but has the advantage of permitting design of the 
filter in both the s and z domains. 

The concept of predistortion applied to the transfer 
function H(s) is introduced. To transform H(s) to H(z) 


the z-trans formation 


z = e ` (5.4) 


is normally used, where 


Ze = sampling time. 
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Because of the cumbersome calculations resulting in the 
use Of (5.4), the simpler bilinear z-transform 


z-l 


_ 2 
Far EM z41 (5.») 





is used. This is an approximation to (5.4) 
The use of bilinear z-transform distorts the frequency 


response characteristics of H(s)| . For this reason 


s=jw 
H(s) is predistorted in such a way that, after the bilinear 
transformation, H(z) will have the desired frequency response 
characteristics. 

Starting in the s-domain a Butterworth fourth order 
band-pass filter is selected because of its flat response 


characteristics. The transfer function of such a filter is 


given by [Ref. 11]: 





bie, "ei og" 
His) = EI ee əəə —— 1570) 
4 o3 b 2 0 4 
s +a o S + (2+ >) ow, cta O s + dÉ 
vhere k = gain 
NM 212, - center circular frequency 


Q = quality factor which is given by Q = w ¿BW and 
BW = 3 dB bandwidth 
a= 1.41421 


(constants given in Tables of Refs. 
5b =r. 0 7 and 11) 
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If BV and Wo are given, the frequencies at the 3 dB 


points, fi and Lor can be defined: 
wg T ој = 212, - 2nf. = BW (5.7) 
ww Np = y (5-9) 


Ihe following restrictions are imposed on the filter 
design using charts. 

(i) The maximum doppler frequency for clutter (angels, 
birds, insects) is 100 Hz. This corresponds to a velocity 
2575107 rum/h for A e .6 m and 22.35 nm/h for A = .23 m. 
The desired clutter suppression at 100 Hz is 30 dB. 

(11) The filter is to be used in air search radar and 
must detect targets with velocities of 881.7 nm/h for 
A = .6 m and 340 nm/h for А = .23 m. This sets the upper 
255 of frequency at 9500 rad/sec or 1572 Hz. 

(iii) Slow moving targets such as helicopters must be 
detected with a suppression of no more than 3 dB. This sets 
the cut-off frequency at this end of the spectrum at 
400 - 600 Hz. 

With these restrictions the following values were 


selected. 


2829.44 4ad/sec, £f, = 445.63 Hz 


dn 1 


15500 rad/sec, f, = 2463.7 Hz. 


E 
| 
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Then, applying basic relationships 


6622.47 rad/sec 


u) 
O 


DIES 2525 


BW 


12616.63 rad/sec. 


with the constraints above and these values the function 
is prewarped and the bilinear z-transform applied. The 
final filter parameters were obtained after successive 
applications. 

In Figure 20 the prevarped form of levee) vs. w 1s 
plotted. The prewarping was made on the center frequency 
wo with the exact Nyquist sampling time, T. - ,2479445x10 ^? 
sec. 

In Figure 21 the final filter power spectrum is given 
after the application of bilinear z-transform. Comparing 
this curve with the analog plot in Figure 19, it may be 
seen that a considerable distortion exists for frequencies 


Ww > Ug. This is reasonable as a result of the application 


of the transformation 


T 
E 2 Ds 
2 = T tan mə (5.9) 
s 
where 5 analog circular frequency, rad/sec 
ün 7 digital circular frequency, rad/sec. 
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For small values of a,, o = Wye but for high values of 
Wr 0.4 ^ Ug. This is implied by the nonlinearity of 


equation (5.9). 


The final filter design parameters are obtained from 


Figure 21 and are as follows: 


SE = 9900 rad/sec 
Е, 2301575703 H7 
ој = 3465 rad/sec 
2| = 551.47 Hz 


Aw = 6435 rad/sec 


Af = 1024.16 Hz. 


More details of the analysis are given in Appendix A. 
Clutter suppression of 34.2 dB was calculated at the 
frequency of 100 Hz. This is better than that of the analog 
delay line which, in the existing MTI radar, is between 


2 nd 30 dB. 
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VI. SUMMARY AND CONCLUSIONS 


This thesis vas motivated by a requirement to find an 
improved anticlutter technique for MTI radars presently used 
in the Greek Navy. Initially, a complete survey of existing 
methods of clutter suppression was undertaken. Table I 
summarizes the results of this survey showing qualitatively 
clutter improvement against various methods and types of 
clutter. 

The Fast Fourier Transform (FFT) appears to be the 
optimum selection. However, it does not lend itself easily 
to introduction into the older radars still in use. There 
is no doubt FFT is used in modern radar systems and will be 
extended in the future, where cost and trained personnel are 
no obstacles. 

For the application under consideration a method of 
digital signal processing was selected. It has the advan- 
tages of being practical, reliable and flexible. In 
addition, due to advances in LSI technology the digital 
filter may be contained in a small package of light weight 
with good performance characteristics and can be constructed 
at reasonable cost. 

The design idea used here was introduced to the author 
for file 71736 timc in a digital filter course at the Naval 


Postgraduate School. 
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The concept of the prevarped transfer function vas 
employed. This approach permits one to carry out his design 
in the familiar s-domain. The design constraints vere taken 
from Table III. This table vas assembled from a number of 
sources giving clutter characteristics for an unusual 
variety of clutter types. Since these -- birds, insects, 
angels and veather (squalls) -- are typical of the ones 
found in the Mediterranean environment against vhich the 
present radars do not give satisfactory suppression. 

The transfer function H(s) vas chosen from a table of 
filters of the Butterworth type. It is a fourth order 
filter. Knowing H(s), predistortion was applied in order 
to compensate for the distortion introduced by the applica- 
tion of the blinear z-transform. This is a powerful digital 
filter design tool since it permits the designer to control 
at each step the values of the filter parameters. 

This thesis actually carried out the preliminary design 
of the selected filter. The results were then plotted 
using Continuous System Modeling Program. The filter 
parameters used in the program were those derived from the 
clutter characteristics surveyed as part of the thesis 
study. At the 100 Hz cut-off frequency this filter gives 
4.2 to 9.2 dB improvement in clutter rejection. This is 
a significant improvement over the existing DLC MTI systems 
and is considered to have met the clutter requirements for 


ships operating in the Mediterranean environment. 
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The adaptive nature of the design proposed for the 
clutter problem suggests that it may be further utilized 
for other applications by interested students in the future. 
Additionally, digital filters may be actually constructed 
using the design program developed here and tested for 


correlation vith the design parameters. 
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COMPUTER PROGRAM 4 


Plot amplitude characteristic of H(z) 
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APPENDIX A 


Analysis of Digital Filter Design 


This appendix covers the analysis in detail of the 
digital filter. The frequency domain transfer function is 


given in Ref. 11, pge 10-4. 





A W AE 
H(s) = O (A.1) 
О) ә 
кәсә то Ыри 4 
О 2 о 
Q 
vhere 


k = 1, gain 


a = 1.41421 
(constants from Table 1, p. 10-4 of Ref. 11) 
ә = p0 
Uu = 6622.47 
(defined 
Q = „525 in the 
analysis) 
BW = 12616.63 


Applying these values to (A.1) 


SS, 1.59x10? s” 
s* 4 1.784x10*s? 4 2.468x10?s^ 4 7.824x10l!gs + 1.923x107 > 
(A.2) 
Applying s > jw (A.2) 
— 27000 0” 70 
Ju; = 8 2 4 3 11 


1277041500 .771.920x107”) 3 3-1. 784x10 u77.823x107-u) 
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This function vas programmed in the computer and the 
results plotted in Figure 19. Associated program vas Program 
Number 2. 

According to the prevarping concept, the center circular 
frequency wo is prevarped. Since BW = 12616.63 rad/sec, 
then Af = 2008 Hz, anā 
3 


_ 1 » m 
T. “Ыл .2479445х10 сес (A-4) 


The distortion factor w' is given by the relation 


Ката (А.5) 





SES E (A.6) 


Applying (A.6) to (A.1) and substituting the known values 


H' (s) = : (u'/Q)* s? : 7 
e 4 1.41421(2)82 4 (24--5) (о) 22 4 5428) 35 + (e!) 
Q Q Q 
(A.7) 
Substituting (A.4) to (A.5) becomes 
w' = 8662.18 (A.8) 
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Since Q = .525, (A.7) becomes 


: E 2.72x109 s” 
H' (s) = — 13 82 12 
s s^ 4 4.2x10 s^ + 1.7х10725 4 5.63x101> 
(A.9) 
Applying s > jw  (A.9) implies 
1. - -2.72x105u” 4 jo 
52 15 кә 12 


070 0 3101:741:7610-“о) 


(A.10) 


This function vas programmed to the computer and results in 
Figure 20. Associated program vas Program Number 3. 
Applying the bilinear z-transform, which is given by 


the relation 





2 z-l 

57 T Z= +Í (A.11) 
s 
to (A.9) becomes 
22 = 22° ae dl 
H(z) = 1 3 7 

35.95869z 45 5.27084z” 4 2.38l4z + 1.03484z + 6.19269 
(A.12) 


Finally, this function was programmed in the compuer and 
the results are plotted in Figure 21. Associated program 


vas Program Number 4. 
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